Orientational correlations in confined DNA 
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We study how the orientational correlations of DNA confined to nanochannels depend on the 
channel diameter D by means of Monte Carlo simulations and a mean-field theory. This theory 
describes DNA conformations in the experimentally relevant regime where the Flory-de Gennes 
theory does not apply. We show how local correlations determine the dependence of the end-to-end 
distance of the DNA molecule upon D. Tapered nanochannels provide the necessary resolution in 
D to study experimentally how the extension of confined DNA molecules depends upon D. Our 
experimental and theoretical results are in qualitative agreement. 
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Thc conformations of biopolymers in living systems are 
often affected by confinement [1, 2]. Examples include 
actin and its analogs in the gel- like cytoplasm [3], DNA 
segregation in bacterial chromosomes [4], and the dense 
DNA packing in eukaryotic chromosomes [5] , in bacterial 
spores [6], and in viral capsids and tail tubes [7-9]. Re- 
strictions of the available conformations fundamentally 
influence function, for example in the case of DNA con- 
densation [10, 11]. 

Single DNA molecules confined in nanofluidic channels 
are a powerful model system for studying the physics 
of confined biopolymers in well-controlled environments 
[12-14]. Understanding the behavior of this model sys- 
tem is thus a first step toward understanding the effects of 
confinement on polymers in more complex biological sys- 
tems. The principal difficulty is that since the persistence 
length of DNA (lp ss 50 nm) is below the diffraction limit 
for visible light, its microscopic configurations are not di- 
rectly observable in the fluorescence microscope. To infer 
the statistics of local conformations it is therefore crucial 
to understand theoretically how such local conformations 
determine large-scale observables, such as, for example, 
the extension R of the confined DNA molecule [Fig. 1(a)]. 

Two generally accepted theories exist for how the ex- 
tension of a DNA molecule varies as a function of chan- 
nel diameter D. De Gennes' scaling theory [15], valid 
for wide channels (D >• £p), predicts that R cx D~ 2 ^ 3 . 
It is based on the notion that a DNA molecule can be 
divided into a sequence of "blobs" [Fig. 1(a)], and that 
the DNA within a blob follows Flory scaling [16]. Odijk's 
theory [17], by contrast, describes the conformations of 
very strongly confined DNA (D -C £p) as almost stiff 
segments deflecting from the channel walls [Fig. 1(b)]. 

Between these asymptotic limits, the local mechanisms 
determining the extension of the molecule are not un- 
derstood, yet this regime is where most experiments are 
conducted. This third regime spans at least an order 
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FIG. 1: Schematic illustration of different regimes for con- 
fined DNA. (a) The de Gennes regime (D > £ P ). (b) The 
Odijk regime (D <§C £p). (c) The Odijk regime with hairpin 
formation (D < lp). (d) A snapshot from Monte Carlo simu- 
lations (described in the text) of a semiflexible chain confined 
to a square channel with diameter D — 4A5£p. 



of magnitude in D, and cannot be understood in terms 
of either asymptotic picture. Several authors have at- 
tempted to understand the dependence of the extension 
upon the channel diameter in this regime in terms of 
power-law relations of the form R cx D~ a [18-20]. In- 
stead we demonstrate that this third regime can be suc- 
cessfully analyzed in terms of the orientational correla- 
tion function C z (si,s 2 ) = {t z (si)t z (s 2 )) ■ Here t z (s) is 
the z-component of the unit tangent vector at contour 
distance s from the beginning of the chain [Fig. 1(b)]. 
Angular brackets denote a time average. A surprisingly 
simple picture emerges: the orientational correlations of 
confined DNA exhibit three distinct behaviors that can 
easily be separated. At short separations the correla- 
tion function decays exponentially. We show that the 
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parameters of this decay depend upon D, and discuss 
the mechanisms behind this dependence. At larger sepa- 
rations, self-avoidance dominates the correlations, giving 
rise to a plateau in C z which can be understood in terms 
of a mean-field theory. This theory applies in the exper- 
imentally relevant range where the DNA conformations 
cannot be described in terms of statistically independent 
blobs [Fig. 1(d)]. Moreover, the theory shows that the 
dependence of the extension upon D is not of power-law 
form for the range of parameters that are easily acces- 
sible in experiments and numerical simulations. Despite 
the fact that the mean-field theory is not based on the 
notion of blobs, it predicts that R oc Z?~ 2 / 3 , consistent 
with de Gennes scaling, in the asymptotic limit of very 
long DNA in wide channels. However, in experiments 
and simulations, the contour length is usually not long 
enough for the asymptotic regime to be reached. We 
show that in this case end effects (quantified by the de- 
cay of C z at very large separations) substantially modify 
the fluctuations of R. We test our model by Monte Carlo 
simulations of a semiflcxible, self-avoiding polymer, and 
by experimentally measuring the dependence of the ex- 
tension R of a confined DNA molecule upon D in tapered 
nanochannels (where the channel diameter D varies grad- 
ually along the channel) [21], that provide the necessary 
D resolution to test the theoretical predictions. We find 
qualitative agreement between the theoretical and exper- 
imental results. 2 A DNA molecule in solution is com- 
monly described as a wormlike chain of persistence length 
£p ~ 50 nm and contour length L, where different pieces 
interact through a screened repulsive electrostatic poten- 
tial. In order to simplify the theory, we instead consider 
a semiflexible chain of N spherical monomers of diame- 
ter a, with a bending potential [^/(fc^T) = — «(tj -tj_i). 
Here = ati points from the center of monomer i to the 
center of monomer i + 1 and tj is the unit tangent vec- 
tor, k is a dimcnsionless measure of the stiffness. Defin- 
ing £p = rea, the wormlike chain can be recovered in the 
limit k — y oc, a — y 0, £p — k& — const. Since the elec- 
trostatic interaction is short-ranged, it can be approxi- 
mated by a hard-core potential with an effective width 
w e g which depends on the ionic strength of the solution. 
Let us constrain the centers of the monomers to a square 
channel of width and height D, extending along the z- 
direction. The end-to-end distance of the chain is given 
by r z — aj^i ^li where if = tj-z (Fig. 1). We characterize 
the orientational statistics of the confined DNA molecule 
by the correlation function C z (i,j) — \tftj). Since the 
confinement breaks the symmetry, the z-component of 
the tangent behaves very differently from the x- and y- 
componcnts. The relevant correlation function is thus 
C z (i,j), as opposed to (t< • tj). The latter function has 
been studied for free DNA [22-24], for strongly confined 
DNA (D < £ P ) [25, 26], and for strongly confined actin 
filaments [27, 28]. 

In our simulations, non- neighboring monomers interact 
by a hard-core potential with effective width w e s = a. 
We use N = 800, re — 8 and a = 1, corresponding to 
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FIG. 2: (Color online) Correlation functions C z (i,j) from 
Monte Carlo simulations, constrained to square channels. 
Left: D = 2£ P = 2na. Right: D = 4A5£p. Regions where 
C z > 0.15 are white. Dashed lines mark the section that is 
shown in Figs. 3(a) and 3(b). 
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FIG. 3: (Color online) (a, b) Sections of the correlation func- 
tion C z (i,j). Blue (upper) line: D = 2lp. Red (lower) 
line: D = 4.45^p. C' z is averaged over all values with 
\i+j — N\ < 10, i.e., between the dashed lines in Fig. 2. 
Dotted lines: estimates of the plateau level C z = (t z ) 2 c from 
Eq. (1). Dashed lines: exponential fits for \i — j\ < n. Dash- 
dotted lines: {tf) c (ij) (hardly distinguishable from solid line 
for D = 2£p). (c) Estimates of the decay parameters a 2 (red 
circles), /3 (blue crosses), as a function of D/£p. (d) (tf) 
[with a(D) and P(D) taken from (c)] as a function of i for 
D = 2£p (blue, upper line) and D = 4A5£p (red, lower line). 
Dotted lines: estimates of {t z ) c from Eq. (1). 

a DNA molecule with w e s ss 6 nm and contour length 
L w 5 |im (assuming £p w 50 nm) . The model and the 
parameters are similar to the ones in Ref. [20]. The 
simulations implement the Metropolis algorithm, with 
crankshaft trial moves [29]. The resulting orientational 
correlation functions C z (i,j) are shown in Fig. 2 for two 
different values of the channel size D, and Fig. 3(a,b) 
shows sections of the correlation functions corresponding 
to the region between the dashed lines in Fig. 2 (left). 

C z (i,j) is seen to exhibit three distinct behaviors we 
now describe in turn. First, when \i — j\ < re, the ef- 
fect of self-avoidance is expected to be negligible com- 



3 



pared to the effect of stiffness. For an unconfincd chain 
(D ^> ^p) the correlation function decays as C z (i,j) = 
a 2 exp{— \i — j\//3} in this region, with a 2 = 1/3, /3 = k. 
The results of our simulations show that the decay re- 
mains approximately exponential for smaller channels 
(D > £ P ), Fig. 3(b). We have fitted a and (3, and find 
that both parameters depend upon D [Fig. 3(c)]. We see 
that a increases as D decreases, reflecting a tendency of 
the segments to align with the channel direction. This 
observation is consistent with recent experimental find- 
ings [14]. The parameter (3 quantifies the initial decay 
of orientational correlations of confined DNA. We expect 
that for D < i-p, the parameter j3 is determined by the 
probability of hairpin formation. The corresponding free 
energy was calculated in Ref. [30]. However, a theory for 
j3 in this regime is lacking, and f3 must be determined by 
simulations. This also applies in the case of the wider 
channels we consider here. 

Second, for \i — j\ ^> /?, the tangent vectors at i and j 
are independently oriented, except for the fact that self- 
avoidance swells the chain. Thus, C z (i,j) — (tf) c (tj) , 
where the subscript defines an average conditional on 
r z > 0. Whereas (tf) = due to the ^-symmetry of 
the problem, (tf ) c takes into account the fact that strong 
confinement breaks the z-symmetry: r z is rarely close to 
zero for a long chain in a thin channel. Figure 3(a) shows 
that this factorization of C' z works very well for thin chan- 
nels. The question is now how (tf) c depends on i. For 
monomers far from the ends, (tf ) is expected to be inde- 
pendent of i (in this case we write (tf) c = (t z ) c ); compare 
Fig. 3(d). As a consequence, a plateau in C z (i,j) devel- 
ops, clearly seen in Figs. 2 and 3(a). We show below how 
the bias (t z ) c (and thus the level of the plateau) can be 
estimated by a mean-field argument. Fig. 2(right) shows 
results for a channel that is so wide that a plateau does 
not clearly develop for a chain of this length. This fact 
is also apparent in Fig. 3(a), solid red line. 

Third, the bias (tf) c is expected to be smaller close 
to the ends, as the chain is more flexible there. This is 
clearly seen in Fig. 3(d) and gives rise to a further decay 
of C z (i,j), seen in Fig. 2 and Fig. 3(a). We note that 
such end decay has been observed in other contexts too 
[22, 24, 26]. However, a quantitative theory for the onset 
and shape of the decay for confined DNA molecules is 
lacking. 

We now discuss the implications of our results for C z 
for the extension of a confined DNA molecule. Let us first 
consider the limit of very long chains (N — > oo). Neither 
our study nor other experimental and simulation studies 
achieve this regime, but it is nevertheless instructive to 
consider. In this limit, fluctuations are negligible, and 
the extension is almost equal to the end-to-end distance, 
which in turn is determined by the correlation function. 
Thus, R 2 « (r 2 ) = a 2 J2i j C z (i, j). Since the contribu- 
tions to the sum from the first and third regions scale 
linearly with N, whereas that from the second region 
grows as TV 2 , the extension is determined by the bias, 
R 2 = N 2 a 2 (t z f c , as N -> oo. We note that in this limit, 



the extension R scales linearly with L — Na, as it must 
when self-avoidance dominates the extension. 

We now demonstrate how the bias (t z ) c (and thus R 
for N — »■ oo) can be estimated by a mean-field argu- 
ment. One way to obtain the correct average (• • • ) c 
for a self-avoiding polymer is by first generating all 
configurations of the corresponding ideal polymer (i.e., 
with spatial overlaps allowed), and then remove all con- 
figurations where two or more monomers overlap. If 
-Pideai(?*z) is the probability distribution of r z for the 
ideal chain, the distribution for the self-avoiding chain 
is given by P(r z ) oc P\& cz \{r z ) A(r z ) , where A(r z ) is the 
fraction of ideal configurations with end-to-end distance 
r z that are free of overlaps. Estimating the functions 
-Pideai^z) and A(r z ) leads to an approximate expres- 
sion for P(r z ). Let us start with Pideal- Unless the 
channel is very thin (D <C ^p), the correlation func- 
tion C^ doal (i,j) decays rapidly to zero for \i — j\ > k. 
Assuming that the correlation decays exponentially as 
Ci doal (i, j) = a 2 exp{-\i-j\/0}, where 1 < (3 < N, the 
distribution Pideai^z) is a Gaussian function, with zero 
mean and variance (r 2 ) « Na 2 a 2 2f3. We now estimate 
A(r z ) by a mean-field argument, similar to the one used 
in deriving the Flory expression for the extension of a 
free self-avoiding chain [15, 16]. If we divide our polymer 
into Np = TV/ (2/3) effective monomers of length 2j3a and 
width a, they are essentially independently oriented. If 
we make the mean-field assumption that these effective 
monomers are uniformly and independently distributed 
within the available volume V = r z D 2 , the probability 
for two given monomers to collide is p = £/V. Here 
£ is the excluded volume of an effective monomer. We 
approximate £ by the value for a stiff rod of the same 
length and width [31]: £ = (7r/2)a 3 (4^ 2 + (7r + 3)/3 + 7r/4). 
Since there are Np(Np — l)/2 w iVj/2 possible collisions, 

the probability of no collisions is A(r z ) = (1 — p)^/ 2 = 
{l-£/(r z D 2 )] N P /2 . With these expressions for P idca i and 
A, we find P(r z ). Differentiation yields the most proba- 
ble end-to-end distance, and thus an estimate of the bias: 

{ )c {^(D)aD 2 ) • [) 

Note that a 2 and (3 (and thus £) depend on D. The 
prediction of Eq. (1) - with a(D) and /3(D) taken from 
Fig. 3(c) - is compared to simulation results in Figs. 3(a) 
and 3(d). The agreement is surprisingly good considering 
the shortcomings of the mean-field theory: First, the ex- 
pression for the excluded volume assumes that the effec- 
tive monomers are randomly oriented stiff rods, while in 
reality they have complicated shapes and have a tendency 
to align with the channel. This assumption overesti- 
mates the bias by an unknown factor in the region where 
a 2 3> 1/3. Second, our theory assumes that monomers 
are uniformly distributed within the channel, whereas in 
fact monomers are more likely to be found in the center 
of the channel than near the walls. This assumption un- 
derestimates the bias by a factor of order unity, for all 
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values of D. 

Equation (1) shows that the dependence of R on D 
is not of power-law form in the limit N — > oo, except 
for very wide channels, where a 2 — > 1/3 and (3 — > K 
[Fig. 3(c)]. In this limit, Eq. (1) gives R oc D~ 2 / 3 , 
in agreement with de Gennes scaling. We emphasize 
that our derivation is not based on the notion that DNA 
within a blob of size D follows Flory scaling, but directly 
takes into account the confinement of the DNA molecule 
in the channel. The theory is valid for channels so thin 
that DNA conformations cannot be described in terms 
of statistically independent blobs [Fig. 1(d)]. Flory scal- 
ing is not relevant to this regime, neither is the question 
which Flory exponent should be used (v = 3/5, 1/2, or 
0.588) [20, 32]. 

For finite values of N the situation is yet more com- 
plicated, as the relative areas of the three regions in the 
correlation function change with D. This is a second rea- 
son why the dependence of r z upon D is in general not 
of power-law form. 

We have verified these conclusions by Monte Carlo sim- 
ulations. Figure 4(a) shows the measured end-to-end- 
distance yj (rf) as a function of D compared with the 
scaling law r z cx D~ 2 I Z . Also shown is (t z ) c obtained 
from Eq. (1), with a(D) and (3(D) taken from Fig. 3(c). 
\J (rf) / L and (t z ) c show qualitatively the same behavior, 
neither of them depending as a power law upon D. How- 
ever, the shape of the curves is clearly different. This is 
a consequence of the fact that the first and third regions 
contribute to the extension - as D increases (at fixed N) , 
the end decay starts earlier, thus decreasing the exten- 
sion faster than predicted by Eq. (1). For even larger 
values of D, the second region vanishes, and the value of 
(rf) approaches that of free chain. 

In order to validate the results from theory and sim- 
ulations we measured the extension of YOYO-labeled A- 
phage DNA (see Ref. [12]). The extension was measured 
as a function of confinement in tapered nanochanncls, 
with a fixed height (D^ = 180 nm) and a gradually in- 
creasing width (from D w = 50 nm to 650 nm). In ta- 
pered nanochanncls the DNA extension can be measured 
continuously as the width increases (for experimental de- 
tails and data for additional ionic strengths we refer the 
reader to the Supplemental Material [33]). For practical 
reasons the experimental conditions could not be fully 
reproduced in the simulations as follows: Whereas the 
simulations were performed in square channels, in the 
experiment the aspect ratio varies with confinement. Ex- 
periment and simulations also differ in that the experi- 
mental contour length is approximately 4 times longer, 
and that the extension R is not identical to the end-to- 
end distance r z . Furthermore, there is no well-established 
consensus in the literature on the effects of intercalators 
on the persistence length of DNA [34-37] nor is the ef- 
fect on the effective width known. Finally, the interaction 
between the DNA and the channel walls could lead to a 
lower effective diameter. For these reasons it is hard to 
compare simulations and experiments quantitatively, and 
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FIG. 4: (Color online) DNA extension from simulations and 
experiments. *J {r'j)/L from simulations (blue diamonds) and 
{t z ) c from Eq. (1) (solid blue line) compared to de Gennes 
scaling (dashed blue line). Extension of A-DNA (48.5 kbp, 
L ~ 20 um [12]) as a function of D = \/DhD w in 180 nm 
deep nanofunnels studied at 0.5x TBE (red crosses). The 
"error bars" show the standard deviation of the distribution 
of extensions, at a given D [33]. Also shown is experimental 
data at 0.5x TBE from Fig. 5 of [18] (black circles). 



we restrict ourselves to noting that they are in qualita- 
tive agreement, as shown in Fig. 4. In particular, the 
extension curve is always steeper than predicted by de 
Gennes scaling, and does not obey a power law. Both 
observations are in accordance with Eq. (1) and are also 
consistent with the observation that the end-decay re- 
gions grow with increasing channel size. 

The results summarized in this paper pose many new 
questions. First, our theory includes the scaling predic- 
tion of de Gennes for D ^> £p , but what does the transi- 
tion to the Odijk regime look like? Second, how do our 
results cross over to the known orientational correlations 
of unconfined DNA [24]? Third, it is necessary to un- 
derstand how the increased fluctuations near the ends of 
the molecule depend on the channel dimension D and the 
properties of the DNA molecule. Fourth, we have seen 
that the parameters characterizing the initial exponential 
decay of the correlation function depend upon D. This 
effect is not quantitatively understood. Fifth, we have 
analyzed a simplified model disregarding possible effects 
of electrostatic interactions between the molecule and the 
walls. These effects are likely to be of importance in the 
experiments we have discussed, and must be investigated. 
Sixth, what does the nonuniform monomer distribution 
of a confined polymer imply for the extension? Last but 
not least it is now experimentally possible to study DNA 
below the diffraction limit of light [38]. In the future, 
these techniques may allow us to directly observe the mi- 
croscopic confirmations of confined DNA. 
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EXPERIMENTAL DETAILS 

Imaging. DNA imaging was made using an epifluorescence video microscopy system consisting of a Nikon Eclipse 
TE2000-U inverted microscope, 60x water immersion (NA 1.0, Nikon) objective (with a 1.5x additional magnification 
lens for a total magnification of 90x) coupled to a back-illuminated EMCCD camera from Photometries, Cascade II 
512. 

Chip design. The measurements were performed in a nanofluidic chip. The imaged part of the chip consists of 
channels with a constant depth of = 180 nm, and a width D w that varies continuously from 50 nm to 650 nm over 
a length of 480 pm. The channels were defined by electron beam and UV lithography followed by CF4/CHF3 based 
reactive ion etching. Finally a thin layer (50 nm) of dry thermal oxide was grown to render the surfaces hydrophilic. 
The devices were sealed with a 500 pm thick borosilicate glass lid. 

Measurement details. All experiments were conducted with A-phage DNA (48.5 kbp, New England Biolabs). 
The DNA was stained with the fluorescent dye YOYO®-l (Invitrogen) at a concentration of one dye molecule 
per ten base pairs. The DNA was dissolved in one of three three different buffers: lx TBE buffer (89 mM 
tris(hydroxymethyl)aminomethane (TRIS), 89 mM sodium borate, 2mM ethylencdiaminctctraacetic acid (EDTA)) 
(Fig. S2), 0.5x TBE buffer (Fig. S3), or 0.05x TBE buffer (Fig. S4). In all cases, 3% ,3-mercaptoethanol was added 
to suppress photoblcaching of YOYO-1. Note that this addition contributes to the ionic strength. Prior to every 
movie recording the DNA was left for 1 min in order to thermally equilibrate after being shifted along the channel by 
applying a pressure gradient along the nanochannel. To control that this equilibration time was sufficient, 5 movies 
were also recorded after 10 min. We did not observe any significant change in the extension of the DNA compared 
to the movies recorded after 1 min of equilibration. Additionally, the relaxation time for A-DNA in nanochannels 
was measured in [1] and found to be on the order of 1 s, which indicates that 1 min of equilibration is enough. Note 
that the gradient in confinement free energy was not sufficient to induce any detectable motion in the DNA during 
the experimentally relevant timescales. 

Data analysis. After equilibration, the molecule was photographed for 400 frames, with a frame rate of 10 Hz. 
From each frame, we extracted the average diameter D — -y/AiAv (defined as the geometric average of the width and 
depth of the channel at the center of the molecule) and the extension r. The resulting extensions from approximately 
100 molecules were sorted by diameter, and collected into bins of 200 measurements each. For each of the bins, the 
average and the standard deviation of the extension was measured. The resulting data is shown in Fig. S1-S4. For fur- 
ther details on the data analysis, and a description of how the extension r was determined for each frame, we refer to [2]. 
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FIG. SI: The average extension r of A-DNA (contour length L « 20 nm) as a function of the average diameter D (geometric 
average of width and depth) in 180 nm deep nanochannels, studied at three different ionic strengths: lx TBE (gray), 0.5x TBE 
(white), and 0.05x TBE (black). 
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FIG. S2: The extension r of A-DNA (contour length L ~ 20 nm) in lx TBE buffer, as a function of the average diameter D 
(geometric average of width and depth), in 180 nm deep nanochannels. The "error bars" denote the standard deviation of the 
measured extension at a given D. 
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FIG. S3: The extension r of A-DNA (contour length L « 20 nm) in 0.5x TBE buffer, as a function of the average diameter D 
(geometric average of width and depth), in 180 nm deep nanochannels. The "error bars" denote the standard deviation of the 
measured extension at a given D. 
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FIG. S4: The extension r of A-DNA (contour length L « 20 nm) in 0.05x TBE buffer, as a function of the average diameter D 
(geometric average of width and depth), in 180 nm deep nanochannels. The "error bars" denote the standard deviation of the 
measured extension at a given D. 



